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Abstract 

The well-known phenomenological dynamic equation for the nuclear shape parameter was 
derived from first principles with the aim of obtaining the microscopic or many-body expressions 
for inertia B, deformation force /, and friction coefficient 7 of collective motion in hot nuclei. 
Nuclear evolution, viewed as a sequence of quasiequilibrium stages, is described with the aid of 
the density matrix pq constrained to given expectation values of nuclear Hamiltonian Ti, number 
operator A/", and operators Q, P, and Ai of coordinate, momentum, and inertia of collective 
motion. Having chosen an explicit expression for Q in terms of the nucleon field operators 
ip^, we construct the corresponding expressions for pq, P, and A4, using a certain canonical 
transformation of '?/'x; "^x' equation of continuity and an assumption that the collective 
variables Pn,t = tr (Pn^tPq) where Pn^t are Heisenberg representations of Pi = Q, P2 = P, 
P3 = Ai, change with time much slower than the number density and the momentum density, 
from which pn,t are built. The same assumption was used to get the closed- form solutions of 
the dynamic equations for p„ j, from which we extract the desired many-body expressions for 
B, /, and 7. After adaptation of those general expressions to a self-confined Fermi gas with the 
phenomenological effective force, they are used for critical analysis of the previous microscopic 
models of those quantities and for elucidating the distinctive features of dissipative collective 
motion in atomic nuclei. 

PACS numbers: 25.70.Lm, 24.10. Cn, 25.70. Jj 

1 Introduction 

Collective or 'macroscopic' motion of atomic nuclei consists in large-scale deformations of nu- 
clear shape, characterized by the parameter q = q (T), which satisfies the equation of motion 

Bq+y%B = f-'jq, (1) 

with q = dq/dT, q = d'^q/dT'^, in which inertia B, deformation force /, and friction 7 are 
certain functions of q. Suggested in [T1121E] for neutron- induced fission, Eq. ([T]), after coupling 
to particle evaporation from strongly deformed nuclei [3] and completing by a random force 
[5l [6], has become a promising tool for an in-depth analysis of experimental data on fusion- 
fission, fusion-evaporation and quasifission processes following nucleus-nucleus collisions at the 
energies per nucleon sizably below the Fermi-energy [3 El El [10] ■ 



The microscopic content of B, /, 7, and q is obscure, because Eq. ([T]) has never been derived, 
as a whole, from quantum-statistical theory of self-bound nonequilibrium Fermi systems. This 
work suggests such a derivation, under hypothesis that Eq. ([1]) is a consequence of the dynamic 
equations for 'macroscopic', 'gross' or collective variables of the system. 

The first step towards microscopic interpretation of ([T]) was made in [HI [12] by postulating 
the microscopic expressions for q and /. The first of those is very simple 



where px is the number density, while 1^ is some function of x, whose integral with the spher- 
ically symmetric distribution equals 0. The / was taken in the form / = —dqW, where is a 
certain functional of px only [12]. Four parameters of this functional: central density, compres- 
sion modulus, central energy density, and a quantity controlling the diffuseness of the nuclear 
surface, were treated phenomenologically. The expressions for W[p] and q[p] were used in the 
variational problem 



whose solution provides the number density p(x. A, p) for the nucleus with the particle number 
A = f (ixp(x, A,p) and deformation q = f (ixy^p(x. A, p). Although Eq. ([3]) closely resembles 
the equation for electronic density in the ground state of a diatom molecule at frozen inter- 
nuclear distance, with W[p] being the analog of the sum of the electronic kinetic energy in 
Thomas-Fermi approximation and the interelectronic interaction, and Al^ being the similitude 
of the electron-nuclei potential, the authors [12] insist that W and py in ([3]) are the total energy 
and number density averaged over a group of levels of a nucleus with fixed A and over A. 

The W in the expression / = —dgW is known as the deformation potential. From the above 
we see that in Ref. [12] it is identified in fact with the internal (or thermal) energy U. This 
contradicts to many other works (see p], and references therein), interpreting W as Helmholtz's 
free energy F either without proof or by assuming that the process of deformation of a hot 
nucleus is similar to compression or decompression of a gas embedded in a heat bath with a 
fixed temperature and occurring so slowly that the temperature of the gas always coincides 
with that of the heat bath. Evidently this analogy is misleading, because the temperature of 
the nucleus in the process of deformation is increasing. 

For cold nuclei the W is identified with the ground state energy of a nucleus with fixed 
A and q. This allows one to obtain it by minimizing the expectation (\E'g| 7i l^^g) of the total 
Hamiltonian Ti. with a state |\E'g) subject to the condition (\E'g| Q |\E'g) = q. Here, Q is the 
operator of collective coordinate 
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where rix = '^/'xV'x is the number density operator and 1^ is the same as in ([2]). To simphfy this 
variational problem, the is replaced by a single Slater determinant |g) of single-'particle' 
wave functions and, to compensate for such a severe restriction of the Hilbert space, the bare 
internucleon interaction in the expression for (gj 7^ \q) is replaced by the effective force. This 
leads to the so-called self-consistent mean-field (SCMF) models. In the early stage [T3], most 
SCMF calculations have used effective forces of the Skyrme type [H] revived by Vautherin and 
Brink [151 [inj- Nowadays [I7j, the three most widely used variants of SCMF's are based on a 
Skyrme energy functional, a Gogny force [18j, and a relativistic mean-field Lagrangian [19j. 

To justify description of many-particle systems with strong interparticle interactions with 
the aid of the SCMF model it is necessary to assume |2Dj that only the two-particle scatterings, 
influenced by the media, are essential, while the three- and more particle collisions are negligible. 
Moreover, the coalescence of few nucleons into a cluster is prohibited. It is hardly possible to 
validate those assumptions for real nuclei. Therefore we trust in the SCMF models simply 
because with several parameters in the effective force, adjusted for a few reference nuclei, they 
reproduce very good the binding energies and densities of many other nuclei |17j . 

Encouraged by good performance of the SCMF theory in studies of these quantities, Flocard 
et al [2T| apply it with a reasonable success to explore the two-humped fission barrier in ^^°Pu 
addressed before within the shell-correction approach [22] ■ Excellent reproduction of fission 
barriers for ^^°Pu and other acitinide nuclei, obtained later with the Skyrme force named SkM* 
[23] , seems mysterious, because the identification of the fission path with a sequence of ground 
states at fixed q is in apparent conflict with the compound nucleus model of fission, used to 
deduce the fission barriers from experimental excitation functions. 

It is well known that for the inertial parameter B in Eq. ([T]) the hydrodynamic formula [3] is 
used by far more often than any other. This poses the problem of the microscopic interpretation 
of this formula, to the solution of which most close is the adiabatic time- dependent Hartree- 
Fock (ATDHF) theory [2HI251I2S]- A nonstationary state vector of the nucleus \qq) = e~*^ |g) is 
generated in this theory by appending to the Slater determinant \q) a phase factor e~*-^, where 
X is a small, time-even, nonlocal one-body operator. Assuming that x has the form x = niqip, 
with 4> = f (ix^x'^x, one expands the energy {qq\ H \qq) in powers of x, stops at second order, 
and dropping the —qm {q\i T-C,(j) \q) term from time reversal symmetry, one gets 

{qq\n\qq) = {q\n\q) + ^Bq\ 

where 

B = m^ {q\i i H,4> ,4) \q) = m j dy. (g| ^x \q) {V4>^ 

has the form of the inertia parameter of an abstract fluid with a mass density m {q\ny^\q) and 
a velocity potential (f)y^. It should be stressed, however, that the expression x = rnqcf) is only a 
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rough estimate of the solution of very involved equations for x- Moreover, the velocity potential 
0x is not proved to satisfy the Laplace equation and it is chosen heuristically. 

Turning to microscopic expressions for friction 7, we remark that for quite a long time 
[271 1211 [291 [Sni [SH [321 [331 [Mj they were extracted from the linear response function of a 
thermally equilibrated Fermi gas confined by the external forces on the small variations of 
those forces. And only recently [35] it has been realized that in the bulk of such a gas, there 
will be no mass flow, jx, defined as the expectation of the momentum density operator 

Px = -(z/2)[CV^x-(VC)^x] (5) 

with the density matrix of the system. Therefore, the collective motion, hence the dissipation 
of collective energy in such a gas will be confined to the area of the nuclear surface. Having 
recognized this fact we have constructed the Rayleigh dissipation function $ from the entropy 
production rate for the canonical ensemble p\ subject to the constraint tr (pxPi) = jx? where 
jx 7^ over the entire nucleus. Although the final expressions for 7 in [35| look reasonable, 
their derivation is far from convincing. First, the p\ was not constrained from the beginning 
by the condition tr {Qp\) = q, instead of which the deformation q was introduced in $ only in 
the final stage of derivation in an intuitive way. Second, having postulated for the mass fiow jx 
the hydrodynamic-like formula: jx = mp^u^, we failed to reveal the microscopic content of the 
'fiuid velocity' Ux, which was treated instead according to the model of Hill and Wheeler [3]. 

The main objective of this paper is the derivation of Eq. ([ID from first principles and 
removing in this way all the obscurities in microscopic interpretations of /, B, and 7, mentioned 
above. Skipping for simplicity the spin and isospin indices and letting h = 1, as done already 
in ([5]), the total energy and number operators are written as 

H = y c^xV-^V^^x + IJ d^dyvil^ - yDr^ryi^yi^., (6) 



Ar = J rfx^:^x, (7) 

in terms of the field operators ip^, ip^, which satisfy the anticommutation relations 

C^y + V'yC = <^(x-y), V^xV^y + ^y^x = 0, CV^; + = 0, (8) 

and a 'bare' NN-interaction f (|x — y|). In order to furnish the formal description of macro- 
scopic motion in a nucleus with total particle number A and total energy E-^, Eqs. ([Sj), ([ZD 
for Ti. and A/" must be supplemented by explicit expressions in terms of tp^, ipy^ for the set 
{Pi, P2, P3) = {Q,P,Ai) of coordinate, momentum, and inertia of this motion. Denoting as 
Pn{T) the expectation values of P„ with yet unknown time-dependent density matrix, we pos- 
tulate that those Pn{T) are very slow functions of T. In formal terms this signifies that the 
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time axis T can be broken up into succession of intervals [Tj, Tj + Tq], where Tq is small enough 
for the expressions 

Pn{Ti + t) - Pn{Td = tdpn{Td/dT, (9) 

to hold for all t C [0, Tq], but yet sufficiently large compared to the internal characteristic times 
to justify description of the nucleus within the corresponding interval of T by the density matrix 
Pq, which maximizes the entropy Sq = — tr (pq In pq) under the following constraints 

tr(Arpq) = A tr(7^pq) = E,, tr(P„pq) =p„(T,). (10) 

According to this definition suggested by Jaynes [36], the pq does not depend on the local 
time t = T — Ti oi the stage in question. The t dependence is attached to the observables: 
Xt = e**'^Xe~'*^, where X are functionals of ip^, ip^- In particular the values of p„ (T) at the 
beginning of the next to i stage are given by p„(Tj+i) = tr (Pn,rqPq), which solves, in principle 
at least, the whole problem of collective or macroscopic motion. 

The operator Q is given in (jlj). The operators P and are expected to be linear functionals 
of the momentum density Px and number density rix, respectively, subject to the conditions 

tr (Ppq) = gtr (A^pq) , tr ([Q, P] pq) = z, 

with q = dq{Ti)/dTi [24J. But those conditions alone do not determine the weighting functions 
in P and in a unique way. Moreover, the equation tr {Ppq) = qti {Aipq) shows that, given 
q, there are correlations between tr {Ppq) and tr (A^pq), whereas the constraining conditions in 
Jaynes's definition of Pq are mutually independent. As shown in section 2, all these difficulties 
in constructing pq are avoided by imposing on pq the following restriction 

Pq = p, [i;'*, , = e-^^'^^x, C = e^^-C 

where the pg is the canonical density matrix, which satisfies the conditions tr {Qpg) = q, 
tr {M Pg) = A and depends on yet unspecified inverted temperature /3, whereas Xx is yet un- 
known real function of the nucleon position x. In order to work out the equations for Xx and 
jS, we use the continuity equation S^px,* = — ""^^^Vjx,* and assume that macroscopic variables 
Pn,t = tr {Pn,tPq) chaugc with time much slower than the number and momentum densities, 
Px,t and jx,t, through which those Pn,t are defined, the argument being that after spatial in- 
tegrations, invoked in those definitions only the slow components Pxj, jxt of px,t, jx,t survive. 
Having postulated explicit expressions for p^^, j^^, we complete the definition of Xx and /5, the 
net result being that the Pq is now completely defined, which allows one to use Eq. ffTOl) for the 
definition of the weighting functions in P and Ai in a unique way. 

In section 3 we utilize another, besides (jH]), an implication of the adiabaticity hypoth- 
esis to get the closed-form solution of the dynamic equation for the expectation tr {Ktpq) 
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of Kt = e'^^^Ke^^^^ , K = P— qAi, and use this solution to express the quotient q" = 
[q (Ti + Tq) — q (Tj)] /rq as a function of the quotient q' = [q (Ti + t^) — q (Tj)] /xq and q = q (Ti). 
After replacement of those quotients with the derivatives, this expression takes the form of 
equation ([1]) for q(T), which permits us to get the many-body expressions for B, f, and 7. 

In section 4, we use the results of Refs. [371 EHl EHl HO] to construct the self-confined Fermi 
gas (SCFG) model, which is a natural extension of the self-consistent mean field model of ground 
state nuclei with fixed shape, to hot nuclei. The expressions for B, /, and 7 obtained for the 
self-confined Fermi gas from our general theory of nuclear collective motion are compared with 
the previous models of those quantities. In section 5 we give a short summary of our main 
conclusions and elucidate the physical nature of collective motion in hot nuclei. 

2 The quasiequilibrium density matrix 
2.1 General expressions for pq 

In this section we define the quasiequilibrium density matrix pq of the system on a separate 
stage of collective motion, which lasts between T = Ti and T = Ti + Tq. This pq does not 
depend on the local time t = T — Ti, while all observables depend on t according to the low 

Xt = e*^*Xe-*^*, (11) 

where X = Xt\^^Q. The values of q(Ti) and dq{Ti)/dTi are denoted as q and q, respectively. 
The total number and the total energy of the system are A and E^. 
Given q, the canonical density matrix of the nucleus reads 

Pg = z;ie-^(«-'^«^+^««), (12) 

= tre-^(«-'^«^+^«'3), (13) 
where the p^ and are functions of P, A, q, determined from the equations 

with respect to p. A, while the inverted temperature /? will be specified later. Putting into (IT^ 
the definitions ([7]), (jl]) of J\f and Q, respectively, and introducing the number density 

Px = Pa,/3(x, q) = p(x; (3, Pg, Ag) = tr (n^pg) , (15) 

we rewrite the equations for pg and A^ in the form 

y rfxp(x; /5, p. A) = A, j dxY^p{yi;f3,p,A) = q, (16) 



from which it follows that 

J d^p^ = A, (17) 
j d^Y^p^ = q. (18) 

The mass flow in pg equals 0: 

tr(p,p,) = 0, (19) 

because Pg, being a function of even operator H — PgAf+AgQ, is itself even under a time-reversal 
transformation, while Px is odd. In order to flnd the general form of pq with tr (pxPq) 7^ 0, we 
use the analogy with the classical fluid with the velocity Ux 7^ 0. According to Mori [H] the 
distribution function of such a fluid, consisting of particles of mass m with coordinates Xj and 
momenta pj, is obtained from that of the fluid with Ux = by the substitutions 

Xi = Xj, Pi p. = Pi - mUx. (20) 

Keeping in mind that fl20|) is the canonical transformation, because x-, p^ obey the same Poisson 
bracket as Xj, pj, and that in quantum systems of identical particles the role of Xj and pi is 
taken over by ■i/'x, "^x; therefore define the density matrix pq with tr (pxPq) 7^ as follows 

Pi = p'g^ P'g = PgW*^^']^ 

where must satisfy the canonical anticommutation relations ([8]) and give rise to the 

primed densities 

< = Px = [cv^:, - (vo<] , (21) 

obeying the conditions 

"-x = "-x, Px = Px - "^^xUx, (22) 

with m being the nucleon mass and Ux an yet unknown c-number vector function. Eqs. fl22l) 
have exactly the same form as the relationships, following from Eqs. fl20|) . between the corre- 
sponding densities in the classical fluid. 

It is readily seen that the required tp'^, ip'* have the form 

C = e-^^-Vx, C = e^^-C (23) 

where Xx is a real function remaining to be determined. Really, the fulflllment of anticommu- 
tation relations and of the flrst condition in fl22l) is evident, while substitution of (123|) into p^ 
defined in ( 1^ leads to the second condition in fl22i) with 

Ux = m- Vxx- (24) 
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The trace of any functional A of the nucleon field operators is 
°° 1 f 

= 5Z M / "^""i' ^^2...c?x;v (0| Vxi^x2-^x^^^:^-CC |0) , (25) 

|0) being the state with no particles. By definition, the primed trace tr'^ is obtained from tr^ 
by putting primes on those ip and which explicitly appear in ( l25l) . but not on those hidden 
in A. From (1251) then follows the formula tr'^ = tr^ , which leads at ^ = Xp^ with accounting 
for pq = p'g to the useful relationship 

tr (Xpq) = tr' (Xp;) . (26) 

As a first application of (!26|) consider the mass flow in p^: 

jx = tr (pxPq) . (27) 

From (l26l) and the expression 

Px = Px + "^"■xUx, (28) 

which follows from fl22]) . one obtains 

= tr' [(p:, + m<u,)p;] . (29) 

By definition, the primed counterpart X' of any operator X, including the trace, is expressed 
in terms of in the same way that the corresponding unprimed operator X is expressed 

in terms of ip:x., "^x- Since, in addition, ip'^, ip'* obey the same anticommutator relations as '?/'x, 
^Z'*, one can drop all primes in (12^ . Then, using (IT^ . (IT^ yields 

jx = mpA,/3(x,g)ux. (30) 

Let us now employ Eq. (I2B]) at X = — In pq to calculate the entropy S'q = — tr (pq In pq) in 
the Pq ensemble. Use of pq = p'g leads to S'q = — tr' (p^ In Pg) , or on removal of all primes, 

S'q = S, 

where S = — tr {pg In pg) is the entropy of the pg ensemble. 

In passing we note that tp'^ and ip'* given in fl23|) can be put into the form 

^'^ = e-'^^ij^e'"', Vx = e-'^^Ce''^, (31) 

where 

X= rfxXxCV'x = / dy^x^n^ (32) 
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will be called the phase operator. This is proven by application of the formula 



e'^Xe-'^ = X + i[A,X] + '- [A, [A, X]] + ... (33) 
at ^ = —X and X = ip^^, and the expressions 

[X, V'x] = -Xx^x, [X, = XxV'x- 

Eq. (|3T1) allows the definition X' = X [ip'*, ip'] of any primed operator to be written as 

X' = e-^'^Xe^^. (34) 

For X = pg on remarking that X' = p'^ = pq, Eq. (IMl) becomes 

Pq = e-'^pge'^. (35) 

This formula for pq is analogous to the definition of the density matrix of the time dependent 
Hartree-Fock formalism [2S1I2S], in which, however, the operator analogous to our x, Eq- fl32p . 
is nonlocal, while that analogous to Pg represents a time-even Slater determinant. 

2.2 The definitions of Xx and (3 

With Eq. (fT2l) for pg, the definition pq = p'g implies that 

= ^-ie-/5(7i'-M.A^'+A.Q')^ (36) 

where Ti.', A/"', and Q' are constructed from tp'^ and tp'* in the same way that Ti, A/", and Q are 
constructed from t/'x and ^Z'*. Use of Eqs. (E]), (JZD, dlD for 7i, J\f, Q, respectively, and Eq. (1211) 
for Ux leads to the expressions 

/in 1 
(ixUxPx + ^ / (^xu^nx, (37) 

M' = U,Q' = Q, so that dSSD becomes 
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Pq = ^ exp I -/3 + A<;(5 - PgA/" - J rfxUxPx + y y c?xu^nx | 



(3J 



where Ux = fn'^'Vxx- In this subsection we formulate the equations for still unknown pa- 
rameters Xx and (3 of pq by exploiting the adiabaticity of collective motion and the continuity 
equation. 

To embody the slowness of collective variables in their formal definition introduce the num- 
ber density px,t at time t 

Px,t = tr (rix.tPq) ■ (39) 
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For t = 0, py^^t reduces to 

Px,o = tr (nxPq) , (40) 

since nx,o = ''^x- Using nx = and the same line of reasoning as in deriving S'q, we find that 
Px,o = (^xPg), or on recalhng ( fTSi) 

Px,o = PA,/3(x,g). (41) 

Having assumed that qt = q{Ti + t) varies with time much slower than px,t, introduce the 
slow component of px,t by the relation p^f = Pa,/3(x, gj. Some collective variables, such as 
qi = tr (QtPq), Mt = tr (A^^pq), which have the following general form J d^Wy^py^^ti where Wx is 
a c-number function, can be put into the form 

c?xpx,iW^x = j dyip^ ^w^ + j dyi5p^^tw^, 

where 5px,t = Px,t ~ Px t is a rapidly changing component of px,t. Now we postulate that 
spatial integration in the second term completely eliminates it, enabling one for estimating 
those variables with the aid of the substitution 

Px,t ^ Px,t = PA/3(x,gt). (42) 

Likewise, the time derivatives of the collective variables of the above type will be evaluated by 
replacing dtp^^t with qtdqpA,i3{y^,qt), where qt = dtqt- 
As seen from the relations 

dtp^,t = dtii (e**^nxe-**^pq) = tr (e^*^^[7^, nx]e"^*^pq) = -m" Vtr (e^*^Pxe-**^Pq) , (43) 
relaying on (15^ . flTT]) and the readily proven expression 

^[H,nx] = -m-Vpx, (44) 

the 9tPx,i can be presented as 

<9tPx,t = -"^~^Vjx,t, (45) 

where 

jx,t = tr (Px.tPq) (46) 
is the mass flow at time t. For t = 0, jx,t reduces to jx given in fl3Up . hence Eq. (HSj) becomes 

^tPx,ilt=o = -V [pA,/3(x, g)ux] . 

If dtPyi,t\i^Q stands under the integral / rfx... over entire space with an appropriate weighting 
function, then we can replace it with qdqPA,/3{'^, q) ■ This yields the equality 

qdqP^ = -V (pxUx) , (47) 
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which we regard as an equation for Ux, because all other ingredients of are known, except 
for P, which will be determined later. Since the left side of ( H71) is linear in q and since 
Ux = m~^Vx:x., (1211), the solution of (H71) has the form 

Ux = gux, (48) 

Gx = V0X, (49) 
where ^x = (qm)'^ Xx is the solution of the equation 

PxA0x + (Vpx) V0X = -5,Px, (50) 

which determines it up to an inessential additive constant. As px is a function of A, j3, q, the 
same is true concerning 0x and Ux: 0x = 0A,/3(x, g), Ux = ua,/3(x, q). 

In addition to using Pxj = PA,^{^,(lt) instead of px,t in the spatial integrals involving px,t, 
the spatial integrals involving jx,t will be calculated with the aid of the replacement 

jx.t -> Lt = WA,/3(x, qt)qtUAA^^ (It)- (51) 



The special form of the slow component jx^ of jx,t in fl5T|) is motivated by the expression 
jx = mp^_^(x, g)gu^^^(x, g), which follows immediately from fISUl) and fHSl) . 

The knowledge of the total energy of the system imposes the following restriction on pq: 
tr (TYpq) = E^. Accounting for fl26l) results in the expression = tr' {Hpg) . Use of (!22l) in 
(1371) leads to the formula 

cixUxPx+— / c/xM^n;^. (52) 

Having inserted this expression for n into E^ = tr' {Hp'^), we remove all primes and invoke 
([H]), (US}, to find 

= - / dxpxM^ + U, (53) 

where 

f/ = tr(7^p,) (54) 

is the internal energy of the nucleus. 
Substituting fHSjl into ([S3D, we obtain 

= MgV2 + t/, (55) 

where 

M^,nL.,^ii. (56) 



Since [/ and M are functions of A, (3, q, Eq. flSSl) is nothing else but an implicit equation for 
P, whose solution uniquely determines /5 as a function of A, E^, q, q. 
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2.3 The definitions of P and M 

Substituting fllSl) into (1551) . we find for pq the final expression 



Pi = exp 

where 



(57) 



Q = j d^Y^n^, (58) 
P = c/xUxPx, (59) 

M=mj dxuln^. (60) 
The expectation values of Q; P; Ai with pq are as follows 

tr (Qpq) = g; tr (Ppq) = Mq; tr (Mp^) = M. (61) 

They are obtained from ([581), iO]), gl]), ([H]); ([59]), ([27D, ([30]), (gH]), ([56]); ([60]), (gOD, dH]), ([56]), 
respectively. 

The operator Q defined in ( 158]) coincides with the collective coordinate Q, introduced in 
([!]). Now we show that operators P and Ai given in (IM?]) and (IBU]) can be interpreted as the 
momentum, canonically conjugated to Q, and the inertia of collective motion. For this purpose, 
we examine the expectation values {Ki)g, {Ki)^ of the operator 

K,=t[Q,P] (62) 

and relate (-P)q with {M.)^- Here and in the following the expectation values of an arbitrary 
observable X in the ensembles Pg and Pq are denoted as 

(X)^ = tr(Xp,), (X)^ = tr(Xpq). 

From (1^ . (155]) . (13^ and the commutation relation 

J dyfyi [uy, px] = -rixV/x, (63) 
valid for any c- number function /x, we deduce 

Ki = i[Q,P] = j cixux j dyYyi [uy, p^] = - J dxu^n^VY^. (64) 
Using (1T5]) yields after a partial integration 

(i^i), = - J cixSxPxVFx = y rfxFxV (pxSx) . (65) 
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The 'continuity equation' (H7|) on accounting for PHI) becomes 

dqP^ = -V (pxUx) , (66) 

so that finally 

{K^)g = -d, j d^Y^p^ = -d,q = -1. (67) 

The second equality follows from f[T5]) . From Eq. flM|) and the identity Uy^ = n'^ it follows that 
Ki = K[, so that, invoking ( 126|) . ( 1671) . we find the equations {Ki)^ = {Ki)^ = —1, which owing 
to (!62|) can be rewritten as 

{[Q,P])^={[Q,P])^ = ^. (68) 

These identities say that in a weak sense the operator P can be regarded as a momentum 
canonically conjugated to the coordinate Q. 

Two last equalities in (1^ lead to the relation (P)^ = {A4)^q between the average values 
of P and Ai at fixed collective velocity q permitting one to anticipate that M = (A^)q will 
represent the microscopic or many-body expression for the collective inertia B. 

While determining Xx and (3 in the previous subsection we postulated that expectation val- 
ues pn,t = Pn(Ti+t) =tr(P„^(Pq) of Heiseuberg representations of collective operators (Pi, P2, P3) = 
{Q, P, Ai) vary with time much slower than px,t =ti'('^x,tPq) and jx,t =tr(px,tPq). This postulate 
was argued by observation that spatial averaging with proper weighting functions performed on 
Px,t and jx,t to get pn,t, remove the rapidly fluctuating components of px,f and jx,t- As a result 
the time t enters Pn^t through qt = q(Ti + 1) and qt = dqt/dt only. Another consequence of the 
hypothesis on slow character of collective motion, expressed by Eq. implies that Pn,t are in 
fact the linear functions of t within the corresponding stage of collective motion. To formulate 
one more implication of the slowness of pn,t, to be employed in finding them, it is convenient 
to put Eq. flFTl) for pq into the form 

Pq = z;'e~^^'^^'>'^^\ (69) 

where 

nc = n + V, (70) 
V = AgQ-qP+^eM. (71) 

The rate of time change of Pn^t is given by dtPn,t = tr (^Pn,tPqj , where P„^f is the Heisenberg 
representation of the fiux Pn = i [H, P„]. Using the cycling property of the trace, the identity 
i [He, Pq] = 0, and Eq. flTTl) . we find 

dtPn,t = tr {i [n, Pn,t] pq) = "tr {i [n, Pq] Pn,t) = tl {i [V, Pq] Pn,t) ■ 

From this expression it follows that the necessary condition for the adiabaticity of collective 
variables is that the quantity i [V, Pq] must be small, which does take place, as seen from 
if i [V, He] is small. 
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3 The dynamic equation for q (T) 
3.1 Collective fluxes 

In this section we intend to derive the dynamic equation for q (T) from first principles. For 
our purposes it is necessary to know the expectations of the fluxes Q, J^, P in both pg and Pq 
ensembles. By definition 

Q = t[n,Q], M = i[n,M], p = t[n,p]. 



Consider first the expectation value ■ Eq. ( 1261) with X = Q in combination with the 
relation Q = Q' — qK[, proven in Appendix A, (11711) . enables us to deduce 



In order to evaluate {^Qj ? we rewrite Eq. ([T^ in the form 

Pg = z;ie-^(^«-^«^), (73) 

where 

Hg = n + AgQ. (74) 

Eq. (1751) shows that TCg commutes with pg. Using this fact and the cyclic property of the trace, 
we find 

(q)^ = tr (t Q] Pg) = tr (t [Hg, Q] Pg) = tr (t [Hg, Qpg]) = 0. (75) 

Substituting (1751) . (IB7|) into (172]) . we finally obtain 



Q)^ = q- (76) 
The expectation value {M.) of A4 = i [H, Ai] is given by 

M)^ = {M)^-q{K,)^, (77) 
where 

K2 = t[M,P]. (78) 

To prove (I77j) we make use of Eq. (|26|) at X = Ai and employ the formula A4 = A4' — qK'2, 
derived in Appendix A, (fT72]) . From ([7S]), dHUD, dSHD, and Eq. dHSD with = ul one finds 

K, = ^[M,P]= ,n j dxdyBj. |„,. pj = -,n / dxn.a.VBj. (79) 
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which on invoking (|T5l) and integrating by parts with accounting for (!66|l yields 

{K2)g = -^j «x5gPx. (80) 

Owing to the commutativity of A4 with Q, we can present Ai = i [H, A4] in the form A4 = 
i [Hg, J^]. Using i [Hg, Pg] = and the cychc property of the trace, we then find that 

= tr (2 [Hg, M] Pg) = tr (Z [Hg, Mpg]) = 0. (81) 

Substituting Eqs. (IHTl) . (IHOl) into (1771) . we finally obtain 

= qm d-xu^dqp^. (82) 



The expectation value ( P ) oi P = i [H, P] can be presented as 



q 

= (^), + ^ (■^), " 2^' ^^'^^ = (^), + 2^'"" y t/x«^9,p.. (83) 

The first equality is obtained using Eq. (!26l) with X = P and Eq. (11731) from Appendix A, 
while the second one follows from (jSTll . ( |80l) . Introducing the operator 

^ = z[7^„P], (84) 

we utilize ( 17^ . ( l62l) to present P = i [H, P] in the form P = — AgKi, which on accounting 
for ( |671) gives 

(^)^ = (•^), + A.- (85) 
Exploiting the identity i [Hg, Pg] = and the cyclic property of the trace, we find 

(^)^ = tr (z [Hg, P] Pg) = tr (z Pp,]) = 0, (86) 

which reduces Eq. (1531) to 

(^)^ = A„ (87) 

hence Eq. (!83|) becomes 

(P)^ = A9 + ^g'm y M^^.p,. (88) 

According to (168!) P is canonically conjugated to Q in a weak sense. Therefore the P = 
i [H, P] can be regarded as operator of the deformation force. Thus, looking at Eq. (18 7p . which 
indicates that the expectation of P with pg equals A^, we may anticipate that A^ will represent 
the microscopic expression for deformation force /. 
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3.2 Microscopic or many-body expressions for B and / 

Now we are prepared to derive the dynamic equation for q{T). The derivation is based not 
on the dynamic equation for Pt = e**'^Pe~**^, as could be expected, but on that for operator 
Kt = e^*^ire-'*^, where 

K = P-qM. (89) 

Since Ai commutes with Q, the K operator is also canonically conjugated to Q in a weak sense, 
but it is easier than P to deal with, because, as follows from the relations tr (Ppq) = Mq, 
tr (Aipq) = M, see (16T|) . its expectation with pq vanishes 

{K)^ = {P)^-q{M)^ = 0. (90) 

The expectation of Kt with pq is designated as kf. 

kt = tr (Ktp^) . (91) 

Differentiating fl9T|l with respect to t, we find 

dtkt = (kt) . (92) 



\ / q 

Using and and remarking that /cq = tr (Kpq) = as seen from (pi . we find that 
kt = (t/Tq)k{Tq), whlch implies that dtkt = k{Tq)/Tq. Substituting this expression into and 
integrating the resulting equation over t with the weighting function e~^*, in which rj = 1/rq is 
the inverted duration time of the quasiequilibrium stage in question, we obtain 

K = l dte-"^' (kt) ^. (93) 

Eq. provides a 'starting pad' for evaluating the dynamic equation for q{T) and obtaining 
the explicit many-body expressions for inertia i?, deformation force /, and friction 7. 

In order to assess the left-hand side of (!93l) . we substitute (l89l) into (!9T|) and use Eqs. (159|) . 
(ISO]) for P, M to find 

kt= I rfxuxj^ , -qm I d^u^^p^^t, (94) 



where px,t = tr (rix,tPq) and jx,t = tr (px,tPq) are defined in fl39|) and (146!) . respectively. In Eq. 
( 1941) both px,f and jx,t appear under the integrals J (ix... over the entire nucleus. Therefore we 
may replace them by their slow components in accordance with Eqs. (1421) . (IHTj) . This gives 

= y c/xUxPa,/3(x, g^Jmg^qUA,/3(x, g^J - gm j dxulpA,f3ix, qrj- (95) 
If we introduce the variations 

5Px = PA,/3(x, g^J - Px, 5g = qr^ - q, = Ua,/3(x, qrj - Ux, 
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then, in the lowest order in those variations, we obtain 

kr = qm / (ixpxUx^Ux + m / dxp^u^Sq. 



Approximating 5ux by (9qUx) Sq, where 6q = qr^ — q, and using Eq. ( l56l) for M , we arrive at 
the desired estimate of the left side of (!93l) : 



^ = j dxp^ (dqul) 5q + MSq. 



(96) 



The following evaluation of the right side of (193|1 is based on a formally exact integral 
equation for Kf. 



Ks,V 



(97) 



which is derived in Appendix B (see also [12] )• By making use of the cycling property of the 
trace and taking into account that Pq commutes with He, we prove for any operator X the 
relationship 

tr (e^*^^Xe"**^^pq) = tr (Xe^^^^^pqe'*^^) = tr (Xpq) . (98) 
Taking the expectations of both sides of flWl) and using fl^JSl) we find that 



Kt) =(K 
q 









K..V 







(99) 



From Eqs. ([89]), ([88]), ([82]) we obtain 



K) =kg-]-q^m I dyMidqP; 



(100) 



so that ( l99l) gets transformed into 



Kt) =Ag--qm I dx.u^dgp^+ j ds [i 



(101) 



Substituting Eqs. ( l96l) . (11 Oil) for k^^, i^Ktj , respectively, into ( l93l) . multiplying both sides 
of the resulting equation hy r] = and introducing the notations 

q" = Sq/Tq = (g^q - q) /t^, q = dq/r^ = (g^^ - q) /xq, 

we obtain 



^qm j dxp^ (dgul) q' + Mq" = Ag - ^q'^m j dxuldgp^ + r(rq). 



with 



roo pt 

r{T^) = ri / dte-"' ds {i 
Jo Jo 



dte 



-r,t 



Kt.V 



J / q 



(102) 
(103) 
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The second expression in fllOSp has been derived from the first one by using the identity rje 
—dte~^^ and integrating by parts. After using the relation 



dgM = m J dx (pxC^gMx + ""x^gPx) 



which follows from ( 156|) . and replacing q' with q and q" with q, Eq. f ll02p reduces to 

qM + ^q'dgM = Ag + r{T^), (104) 

where M and A^ are defined in ( 156|) and ( fT6|) . respectively. 

The composition of Eq. (11041) would be fully identical to that of the dynamic equation ([1]) 
for q{T), if r(rq) in (11041) were linear in q. In subsection 3.3 we show that in the adiabatic 
approximation this is the case. Therefore we find the following many-body expressions for 
inertia B and deformation force / in ([T]): 

B = M, f = A„ (105) 

which confirm our expectations concerning M and A^. 

In order to find an alternative expression for /, we prove the formula 

Ag = -dgF, (106) 

where the F is Helmholtz's free energy: 

F = U- (107) 

As a first step we use Eq. f|T2|) . the normalization condition trp^ = 1, Eq. fl54|) . and equations 
f|T4l) at yU = fig, A = Ag, to put the entropy S = — tr(pg In pg) into the form 



S = \nZg + (3U - pfigA + pAgq. (108) 

Next we put this into (I107P and differentiate the result over q, which yields 

- d,F = Ag + p-'dg In Zg + qd.Ag - AdgPg. (109) 

From Eq. f|T3l) and equations in f|T4l) with fi = fig, A = Ag we derive the relation 

f3-'dg\nZg = /5-iZ;Hr9,e-^(^+^««-^''^) = -tr [pgdg (A,Q - pgAf)] = -qd.Ag + Ad.pg, 

which reduces ffT09|) to ffT06|) . 
According to the expressions 

/ = -d,W, f = Ag, Ag = -d,F, 

the deformation potential W must be identified (up to an additive constant) with Helmholtz's 
free energy F. This justifies the practical procedures of calculating W in hot nuclei described 
in [B], and the references cited therein. 
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3.3 Many-body expressions for 7 

In order to facilitate the analysis of r(rq), the flux Kt in fll03p is rewritten in the form 



with 



(110) 



:iiii 



= I [n„ K] . 

Eq. ffnol) follows from Eqs. ([TOD, ([89]), [Q,M] = 0, ([62]), ([78]). Since pq commutes with He, we 
have {n = 1, 2) 

tr (^ [K, K^,t] Pq) = tr (z [He, K^,tPq\) = 0. (112) 

The last equality follows from the general relation tr ([X, F]) = due to the cyclic property of 
the trace. Combining fITO]) with (11121) allows to write: 

tr {t [Kn,t, V] pq) = tr {t [Kn,t, K - H] Pq) = tr (K„,tPq) , (113) 

where Kn,t = i [H,Kn,t]- Substituting (II 101) into (11031) and accounting for (11131) . we obtain 

1 



An 



(114) 



According to Appendix C, the integrals dte~^^ {kn^tj ,"^ = 1,2, are proportional to g, 
which is a small quantity of our problem. Therefore the second and third terms in (11141) lead 
to the corrections to A^ and |g^(?gM in (I104p . which can be neglected. With this in mind we 
find 



dte'^'Hi 



Pq 



or, using k^ = e'^^^k^e and the cyclic property of the trace 



c/te-^r e'*^K'=e"'*^i [Pq, ^] ) • 



;il5) 



Regarding as a small perturbation of Tic in the expression 7i = Tic — V enables one to replace 
the operators e^**^ = e^***^''^''"^'' in (IllSp with e^**'^% because an error introduced into r(rq) by 
this replacement is of second order in V. This gives 



dte-'^hi e'*^^ir"e-**^^i [pq, V] 



dte-"^' 



■ (116) 



The second expression is obtained from the first by cycling permutation of operators inside the 
trace. According to Eqs. (11741) . (11751) . (I176P of Appendix A, the He, V, k^ are related to the 
primed operators as follows 



He = H' V = AgQ' - qP' - -rM', K'' = T'. 



(117) 
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Inserting these expressions into (11161) . accounting for (|26|1 . and removing all primes, we find 

1 



r(rq) 



dte 



T , e 



-UHn 



AgQ -qP~ -q'M e 



JtHn 



118) 



According to flllSp r(rq) contains in addition to the term proportional to q also the terms 
proportional to and g^. Evidently, these additional terms, upon putting (11181) into (I104p . 
will modify A^ and ^q'^dgM. Since these modifications would destroy our expressions (11051) for 
B and /, whose validity is beyond doubt, we require that the coefficients in front of the terms 
proportional to A„ and q^ are negligible small: 



<1, 



POO 

Jo ' 



< \dgM\ . (119) 



As a result, we obtain the dynamic equation for q{T) in the final form: 

1 



where 
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qM + ^q'^dqM = Ag - qi, 



(120) 



(121) 



The physical condition for the validity of (I119P is the adiabaticity of collective motion. To 
show this we note that adiabaticity implies that Q and Ai are 'quasiintegrals of motion', which 
means that their commutators with Ti, which are equal to those with Tig, are small. Letting 
[ng,Q] = [ng,M] = 0, we find that e-^*^«ge^*^« = Q, e''^'^^ Me'^^^ = Min the static limit. 
The identities ([JF, Q])^ = 0, ([JF, M.])g = 0, proven in Appendix D, then show that in this limit 

dte-''^' {[T,e-''^^Qe''^^]) =0, / rfte""* ( [j^, e-^*^«A^e^*^«] ) =0. 

'0 ^ Jo ^ 

These identities confirm that for sufficiently slow collective motion the conditions (I119P are 
satisfied. 

Inserting into (11210 the identity 1 = dtt, integrating by parts and using the relations 



-itn 



itn„ 



-itHn T^itHn 



the last of which follows from (IHH), we obtain 



7 = ^ 



noo poo 

/ dte-''H{i[T,e-''^^Pe'''^^]) + / rftte"^"* (i [j^, e-'*^«J^e'*^«] ) . (122) 
Jo ^ Jo ^ 



In the classical case one describes the nucleus in terms of the coordinates Xj and momenta pj 
of the particles j = 1, ...A, bound by the mean field potential V{yi, q). The P and T = i \Hg, P] 
are given in this case by YlVj ^^'^ Yl respectively, with pj = Ux^Pj, Fj = — Ux^. Vy(xj, q). 
Recalling the definition rj = I/t^, we see then that the first term in (I122p is smaller than the 
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second one by a factor of T/xq, with T ~ YliPjl Yli^j denoting the characteristic period of 
particle motion in the direction of the mass flow. As far as characteristic times of the internal 
motion are much shorter than Tq, we may omit the first term in (11221) to obtain 

/»oo poo 

Jo Jo 
Now we use in the second term the cyclic property of the trace, the commutativity of Tig with 
Pg and the substitution t ^ — t to join it to the first one: 

/oo 
dte-'^l*ltC^^(t), (123) 
-oo 

where Cy^-f (t) is the time correlation function of JF: 

(t) = tr (e^*^«^e-^*^«^pg) = tr (e**^«^e~**^«^Pg) . (124) 

Note that the Hamiltonian Tig, generating the time evolution of the first in ( 1124p . dif- 
fers from the true Hamiltonian by a term AgQ, which is the analog of the electron-nuclei 
interaction component of electronic Hamiltonian in the diatom visualization of nuclear collec- 
tive motion. The Hamiltonian Tig enters Cjrjr{t) not only through pg and e^**^», but also 
through T = i [7ig,P], in which the P itself can be put into the form P = ■mi[Hg,(f)\, where 
4> = J dx(f)xn:xi- This expression for P is proven by using Eq. (H9l) for Ux in Eq. (1591) for P, 
integrating by parts, and employing the relation Vpx = ~mi[Hg,ny_], which follows from Eq. 
(jH]), because of commutativity of Q with n^. Using the above expressions for JF and P, we 
put Cjrjr (t) into the form 

Ctt (t) = m'^^, % it) = tr (e^*«^</>e-^*^«0p,) . (125) 
Combining Eqs. fll23p . fll25p . we obtain the expression 



7 = irri^ 



/oo d o 1 1 ) 
dte-'^\hC'^ (t) , C'^i (t) = (126) 



relating friction with the time correlation function of the reduced phase operator (p = [qm) ^ x 
in the pg ensemble. 



4 Collective dynamics of self-confined Fermi gas 
4.1 Deformation force 

For practical calculations of key ingredients (/, B, 7) of equation ([T]) we treat the nucleus 
as a self-confined Fermi gas of quasiparticles bound by a self-consistent potential generated 
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by the effective force. The quasiparticle energies are denoted as Ex, their real- valued wave 
functions being ipx (x). According to [371 EH], the entropy S and the one-body density matrix 
Pyz = {i^yi'z) in this gas have the same expressions in terms of occupation numbers fx as for 
noninteracting fermions 

^ = - E 1^ ^ + (1 - ^) 1^ (1 - ^)] ' (127) 

A 

Pyz = $^/AV^A(y)¥^A(z), (128) 
A 

where 

= (l + e^^^^-'^))"'. (129) 

The Gibbs grand potential Q [p] of the system with Hamiltonian Tig and density matrix p, is 
given by 

n [p] = tr [{Hg - pAf) p] - p-'S [p] , (130) 

where S [p] = — tr(plnp) is the entropy of the p ensemble. A remarkable property of ^[p], 
established by Gibbs in the classical case, is that it reaches its minimal value flmm for P = Pg- 
For quantum case this property is proven in [39] and [IQ]. It follows from fllOSp . fl74|) . fl5^ . f[T^ 
that 

^min = [Pg] = "/^"^ In Zg. 

Let p differs from pg by small variations of wave functions ipx (x), but has the same values 
of P, p, A, Ex as in pg. Since /J'^S* [p] for so defined p coincides with l3~^S, because S accord- 
ing to fll27p . fll29p is not depending on (px (x), the variational derivative of Q [p] over ip^, (x) 
coincides with that of tr [{Hg — pAf) p] . Using Eqs. fl71|) and ([7j) for Hg and Af, respectively, 
the variational equations for (p^ (x) can consequently be written as 
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U[p]+aJ rfyYyPy -pj dypy - ^ fxex j dy^l (y) j = 0, (131) 



S(Pu (x) 

where U [p] = tr (Hp) and 

Py = E^^A(y) (132) 

A 

are the internal energy and the number density in the p ensemble, respectively, while fxCx are 
the Lagrange parameters to ensure that each wave function (fx (y) should remain normalized 
under variations. 

In order to present U [p] = tr {Hp) as a functional of ipx, we first use Eq. ([6]) for H to write 
U [p] = [dy (ry^i^y) + \l dydzv{\y - z|) (V^^CV^z^y) , (133) 
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where (...) = tr (...p). The effective interaction ^'''''(y, z) is defined by the relation 

j dydzv{\y - z|) (z^*V^*V^2Vy> = j dydzv^^{y, z) (pyp^ - PyzPzy) • (134) 
Then on using ( ]128p . Eq. ( 1133^ becomes 

^ = i E / ^y^^ [^^^ (y)]' + 1 f dyd7.v''''{y, z) (pypz - PyzPzy) . (135) 

Assuming that f'^^(x, y) does not depend on pn, where R =i(x + y), we find from (I13ip . 
(11351) the infinite system of equations 

- T^VV. (x) + + AY^) (x) = E,^, (x) , (136) 
2m 

where = + p. and 

l^HF^^ (x) = j dy^;<=ff (x, y) [py^, (x) - p^y^, (y)] , (137) 

which constitutes, in combination with equations (11281) . (11321) for py and Pxy, the formal basis 
of the temperature-dependent Hartree-Fock method. 

In order to account for the fact that f'^^(x, y) does depend on pR,, it is necessary to specify 
the exact form for such a dependence. For instance, in the case, when f''''^(x, y) is the Skyrme 
force t>sky(x, y), Eq. (I135p takes the form Usky = J c^xi^sky (x), in which -f^sky (x) is a certain 
algebraic function of the number density px, the kinetic energy density 2mrx, the spin-orbit 
density Jx, and their derivatives [23j. The equations for ip^, (x) will then be the same as ( I136p . 
while ( I137P takes the form 

ViS^x (x) = V/i V(/.A (x) + f2^x (x) , (138) 

where /i depends on px only, whereas /2 is an algebraic function of px, Tx, Jx, and their 
derivatives. 

The solutions fu of Eqs. (I136p . (11370 with py, pxy defined in (I132p . (I128P depend on 
P, p. A: Px = p(x;/3,p. A), pxy = p (x, y; /?, p. A). Given px = p(x;/?, p,A) the chemical 
potential fig and deformation force / = A^ as functions of f3, A, q, are determined from Eq. 
(1161) . Expansions for and Jx in powers of h, involving px and its spatial derivatives, which 
are converging at arbitrary x, may be obtained at finite temperatures only [23l |13] . Using few 
lowest in h terms of those expansions in i^sky (x), one gets the so-called extended Thomas- Fermi 
(ETF) internal energy density H^^^ (x), which is a function of px and its spatial derivatives: 
-^sky^ (x) = -f^sky^ (Px^Vpx, .•). The internal energy in this case becomes a functional of px only, 
which allows to use the extremum property of Q [p] to write the variational equation for px 

6 



5p^ 



j dyHi^/ (Py,Vpy, ..) + A j dyVyPy ~ fx j dypy^ = 0. 



(139) 
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This equation has the same form as the Strutinsky-Tyapin equation ([3]) commented in Introduc- 
tion. According to [23], Eq. f ll39p enables one to calculate px and U two orders of magnitude 
faster than the HF formalism. Moreover, the potential = V/iV + /2 obtained with ETF 
values of px, ^x, Jx, being used in the HF equations, provides a rather inexpensive and yet 
accurate shortcut to fully quantum HF calculations. 

For P"^ — > +0, the f\ goes to the step function and Eqs. (11361) . fll37p convert into the 
system of A usual HF equations, providing a formal basis of self-consistent mean field models. 

As stressed in [H] and [121 US], in contrast to the liquid drop model, the HF formalism 
and its semiclassical version, Eq. fll39p . are capable to account for though small but finite 
compressibility of the nucleus. In order to illustrate this, consider the density distributions px 
obtained in ETF calculations [1^ W7\ performed for ^^°Yb with a realistic Skyrme interaction, 
namely SkM* [23]. The function was taken in the form = Iz^ — r^, with r = \/ x\ + 
z = X3, where Xj are the Cartesian coordinates of x. The collective coordinate q for such a 
choice of Yy^ has the meaning of the quadrupole moment q = Qi- Since Fx is invariant under 
rotations of x about the Xj, = z axis, the densities Px are functions of 2, r: px = p{z,r). The 
values of p(z, r)/po, with po = A/{^7rRl), Rq = 1.18 A^/^ fm, along the fission path of ^^^Yb 
are depicted in Fig. 1. The profile function Pc{z) is defined as such contour on the z, r plot 
of p{z,r), where p(z,r)/po = 0.5. In Fig. 2 one can see those Pc{z) for the sequence of p(z, r), 
presented in Fig. 1. Let V = V {Q2) be the volume of a body of constant density po bound 
by the surface of revolution of Pc{z) about axis z. From Fig. 3 one can see that this volume 
slightly decreases with elongation, which means that on its way to scission the nucleus is getting 
denser. Using the empirical formula (see Eq. (2.70) in [48J) 

R ^ (1.12A^/3 - 0.86^-^/3) fm 

for nuclear radii R, we find that the ratio of the total volume 2Vf[ of two spherical fission 
fragments to the volume Vcn of the spherical compound nucleus ^^°Yb is 2Vs/Vcn = 0.954. This 
number is surprisingly close to the value \^(240b)/^(0) ~ 0.95 from Fig. 3. 

4.2 Inertia 

By using (I132p . Eq. (!56l) for inertia M may be rewritten as 

M = mJ2fx I rfx<^Ux)M^ = m5^/,(A|M^|A), (140) 

where is the one-body operator with the matrix elements (x| |y) = m^5(x — y). Since 
Ux = V0X, the inertia is expected to be expressible in terms of (pxu = / d^Vx (x) ^xV'i^ (x). The 
derivation of such an expression begins with using (11290 . (I136P to represent (11400 in the form 

M = m^(A|u2(l + e^(^-^))"V)- (141) 

A 
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Here, H is the single-particle Hamiltonian 



with p = —iW and 

V^ = V^^+KY^, (142) 

in accordance with (11361) . (I138p . In the subsequent discussion the V/iV term in is for 
simplicity omitted. Now, using the identity 



we obtain 



and 



—I 
2m 



-I (V0x) , 



V0X P + P V0: 



[[H 



1 Y'xJ ) Y'xJ 



If ~ X2 

m V 



Next, replacing V0x in (I145p with Ux in accordance with fHUl) . we have 



1 S^xj) S^xj 



(143) 

(144) 
(145) 

(146) 



Substituting f ll46p into f ll4ip and accounting for f ll36p . we find the desired expression for M: 



(147) 



where Exu = Ex - E^. 

For calculating M with the aid of (11470 one should solve Eq. (!50|) for 0x) which in the nuclei 
of complex shape is quite tedious. The well-known simplification consists in using the Werner- 
Wheeler approximation for Ux. This approximation does not satisfy the condition Ux = V0x, 
hence there is no Werner- Wheeler version of the phase function 0x, to be used in (I147p . To 
circumvent this obstacle we express M via the matrix elements of the one-body operator 



F = i 



H, ^(uxP + pu. 



:i48) 



Replacing V0x in the right-hand side of (I144p with Ux, we find the formula 



-(UxP + pUxj 



which allows to put Eq. (11480 for F into the form 



-mlH, 



F = —m 



H, [H, 
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from which, using (11361) . we find 

Fx, = -mElJxu- (149) 
Inserting (11491) into Eq. fll47p . we arrive at the required expression for M in terms of Fxu'- 

M = -J2\Fxuf{fx-fu)/El. (150) 

Xu 

In order to use f llSOp in practical calculations, Eq. (11481) for F is rewritten in the form 

p = piv)^p{u)^ (151) 

with 

F(^) = -GxVKx (152) 

and 

where u^^ = d {ux)i /dxk- The Latin indices represent the Cartesian components of the corre- 
sponding vector. The repeated indices are summed over from 1 to 3. 

Following p9], consider a simplistic situation, when Ux and are given by: 

1) Ux = if X is off the nuclear surface and Ua = qug, if x = a, where a is an arbitrary point 
of the nuclear surface and na is a unit vector normal to the surface at point a; 

2) Vx. = —Vq safely inside the nucleus, = safely outside of it, and is a. step-like 
function of ^na with vanishing width for x close to the surface, where ( = (x— a) na and a is 
that point on the nuclear surface, which is the most close to x. 

It can be shown, using Eqs. ( HHl) . (11511) . (11521) . (I153p . that F goes to dqV^ in this case, so 
that M given in (I150p goes to 

M- = -Y, \id,V)xuf {fx - U) /El. (154) 

The index cr is designed to stress that Eq. (I154p . according to [31], follows from the cranking 
model formula for inertia of a Fermi gas [27] , if any residual interaction is neglected. It can be 
shown analytically and numerically calculations [21] that M^"^ as a function of q is divergent 
near the crossing points of the single-particle energy levels, where Exu = Ex — = Q. This 
signals that the simplistic situation considered in ^29j is physically meaningless. Our derivation 
of (11540 shows that M'^^ could have been applicable for such a hypothetical nucleus, in which 
the mass flow induced by the distortion of the nuclear surface is getting immediately damped. 
This implies, of course, a very short mean free path, which is incompatible with the Fermi gas 
model of the nucleus used to derive (I154p . 

The derivation of M^^ given in [23 EI] is based on the linear response theory. The fact 
that this theory leads to wrong inertia means that its application to nuclear collective motion 
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is inappropriate. Really, the LRT regards such motion as a sequence of the states of complete 
thermal equilibrium at fixed nuclear shape. Within this simplistic picture there is no room for 
the mass flow in a hot time-evolving nucleus, which precludes the possibility for LRT to get 
the realistic estimates not only of collective inertia, but also of friction (see [35] for details). 

In passing we note that some formally erroneous simplifications, destroying the self-consistency, 
need also to be made in deriving the cranking formula in cold nuclei from the adiabatic time- 
dependent Hartree-Fock theory (see Section 9.3 of [25] for details). Though the cranking model 
expression for the inertia in the dynamic equation for q is wrong, the cranking model formulae 
for the inertia of the rotational and translational motions are correct. This is because the op- 
erators, which generate rotations and translations commute with the total Hamiltonian, while 
the operator P of collective moment conjugated to Q does not. 

4.3 Friction 

Let a^, ax represent the creation and annihilation operators of quasiparticles with energies Ex 
and wave functions ip^, (x), as found from Eqs. (11361) . (11381) . Substituting the expressions 

ilJ^ = ^a^ip^{yi) , = ^ a^(^A (x) (155) 

into (11261) and using (11251) . we find 

^ = irr? ^ 4>\u4>x^u^ / (ite"''l*lt— (a^(t)a^(t)a^^a^^) (156) 

where 

For practical calculations of (a^ (t) a^ (t) a^^aj^^) we take Pg in the form 



Pa = ^/ exp 



-P ^ (^A - /i) a*xax 



The Hamiltonian Tig, which enters the propagators in C^^ (t), is modelled as 

Hg = Ho + 

where 

Ho = Y^Exalax = j rfxT/^^V^x + j d^r^V^^;^ (157) 

A 

is the mean-field Hamiltonian, while V^'^^ is the residual interaction. The second equality in 
(11571) is obtained by using (I136p . (I155p . The V^^^ is given by 



^^a*a 
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where 



^A^^ = j dyidyifx (x) f '^'^(x, y) [pyip^ (x) - p^yip^ (y)] . 
Using ffT32|l . ffT28D . we obtain 

=^f» rfxrfy^A (x) (y) ^^''^(x, y) [(p^ (x) (y) - (x) (y)] . 

Let to represent the colhsion time: to ~ tq/vf, where ro is the range of the effective force 
and vp is the Fermi velocity. Then, following ^91131], we postulate that for |t| > to 



(«A (^) «^ W «Ai«<^i) = («A (^) (^) all) + ^^a/a5i.iAi/ai, (158) 

where 

(a* it) a,,) = 5.,A/Ae~^^^l*'e^^^*, (a, {t) a^) = 5.a J.e-^^^l^'e"^^^*, (159) 

with fy = l — fy. Within perturbative treatment of V^''*''^, the collision rate Tx of quasiparticles 
with nuclear media is related to the effective force as 

Ti = 27r ^ (12| v^"" |34) (43| v^^ |21)^ 5 (E1234) {f2hh + hhU) , (160) 

234 

where 

= 1-^1) ~ ' -^1234 = El + E2 — E3 — E4. 

The bold digits for quasi-particle quantum numbers are used to stress that the effective inter- 
action v"^^ in Fi is not assumed to be diagonal over the spin-isospin quantum numbers s. If v"^^ 
does have this property, 

{SiS2\v''^{x,y) \S3S4) = 5siS3'^S2^4^(x,y), 

and then if w(x, y) has the form 

f(x,y) = Vo(x)(5(x-y), 

then Fi becomes 



Fi = 67r ^234(5 (^1234) (/2/3/4 + /2/3/4) 



234 



"^1234 = j dyiipi (x) ip2 (x) Vo(x)v93 (x) (fi (x) . 

The factor Gvr instead of the usual 2tt arises due to the fact that for the s independent zero-range 
interaction the summation over S2, S3, S4 is reduced to finding the sum 

(S1S2I IS3S4) (S4S3I (1 - p^'^) \s2s1) = 3, 
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where P*-*-* is the spin-isospin exchange operator. 

Eqs. (11581) . (I159P imply that for |t| > to the trajectory of any selected particle does not 
correlate with the trajectory of any selected hole, unless their quantum numbers coincide. 
Both the particle and the hole, however, have a chance to collide with nuclear media, which is 
reflected in the damping factors in Eq. (1159^ for their propagators. The simplest derivation of 
Eqs. (I159p . (11601) . not involving a diagrammatic expansion, but using instead Tserkovnikov's 
substitution for {ci*x (t) a^-^) ^ and (^ay{t)a\^^, can be found in [33]. In [M] we use a simple 
version of the effective force to obtain the semiclassical expression for Fi in finite spherical 
nuclei and compare it to the values of Fi in infinite matter. 

The smallness of the collision time to compared to the mean free path time ~ F~^, allows 
to use the asymptotic form (11581) . (11591) of (^a\ (t) (t) a*Xj(^vi)g for all t in (11561) . giving 

where E^x = Ey — Ex, Txv = \ (Fa + F^,). Then, we use the formula 

/CO t4 



to obtain 



^ = -4m'J2<PluExurxJxfu 



Noting that the (pl^^Exu^xu factor in the summand is antisymmetric under interchange of A with 
z/, we replace the remaining factor fxfu by its antisymmetric component 

{fxfu)^ = 2 {fxfu - fufx) = 2 [/a (1 ~ fu) - fu{^ - fx)] = 2 (/a - fu) 

to obtain 

7 = -2m'T J2 (/a - fu) Exu, (161) 

uX 

where F is an 'effective' value of F^^a, not depending on the quantum numbers. 

A seemingly different expression for 7 follows from (I123p . if we replace there !F with J^q = 
YliXu ExuO'x'^ui where Fxv = (A| F |z/), with F defined in (I148p . The replacement of JF with JFq 
is motivated by the assumption that V'^'^^ is small and it leads to the quantity 

2Ex,T 



7o 



uX 



(El + F2) 



From Eqs. (I15ip . (I152p . ( I153p for F it is seen that Eq. (11620 coincides with the formula of Ref. 
[35] . Since the assumption on the smallness of V^'^^ implies that F is small too, we may omit F 
in the denominator of (11621) . which on accounting for Eq. (I149p for Fxu leads to an expression 
for 7o, which is identical to (11611) . 
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Comparing Eqs. (I16ip and (11471) for 7 and M, respectively, we find the formula 

7 = 2TM, (163) 

which shows that the q dependence of the reduced friction coefficient 7/M is completely deter- 
mined by that of the collision rate F defined in (I160p . For practical use of this expression we 
rewrite it as 

7 = 2\~^vfM, (164) 

where X = Vp/T is the mean free path of the nucleon. 

In Figs. 4, 5 the 7 from (I164p are superimposed on the semiclassical values 79'^ of 70 obtained 
from (11621) in [35j by a conventional Monte-Carlo method. The calculation is done for ^osp]^ 
and ^^^Ds at A = 20, 50 fm along the fission path starting from the spherical shape with the 
radius R = 1.12^4^/^ fm and ending at a necked-in configuration with the neck radius of 0.3-R. 
The M given in fl56|) was calculated for px = 0.17 fm~^. For q we took the distance between 
the mass centers of the two halves of the nucleus. The Werner- Wheeler expression for the 'fiuid 
velocity' Ux and the Cassini oval parametrization for the profile function Pc{z) of the nuclear 
shape were used. One can see that 79"^ does reproduce the general trends of 7 as a function 
of A, A, and q, although there are differences in detail. This may be related to the fact that 
the summand for 7q'^ is a sign-alternating function of the initial position and velocity of the 
particle, which makes the direct trajectory-integral simulations of 7q'^ unreliable (e.g. see [50]). 



5 Summary and some remarks 

The main objective of this study consists in working out a microscopic description of collective 
motion in hot nuclei, allowing to express the parameters of a phenomenological collective model 
[H [21 E] in terms of nucleonic quantities. The study starts with the introduction of the oper- 
ator Q of collective coordinate and solving the two interwoven problems: construction of the 
quasiequilibrium density matrix pq and establishing the explicit expressions for the operators 
P and Ai of collective momentum conjugated to Q and collective inertia, respectively. 

The crucial steps are: 1) the substitution t/'x e~*^''^/'x5 "^x ~^ e^^^^ip^ for the nucleon 
field operators made to get pq from the canonical density matrix pg subject to the condition 
tr (Qpg) = q and 2) the conversion of the adiabaticity hypothesis into assumption that the 
deviations of the number density Px,i and the momentum density jx,t, which are rapid functions 
of time, from their slow components Pxt and jx^ defined in ( H2l) and ( |5T1) . respectively, are 
smoothed out to a negligible level after performing the x integrations invoked in the definitions 
of collective variables tr (QtPq), tr (PtPq), tr (A^^pq) in terms of px,j and jx,t. 

Equating two different expressions for the time derivative of tr {KtPq) with Kf = Pt — qM.t 
then led us to a relation between the variation 5q of collective velocity over the quasiequilibrium 
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stage in question and a time integral over this stage of the function tr (z [H,Kt]). After esti- 
mation of this integral in the adiabatic approximation, the relation takes the same form as the 
equation of motion of the phenomenological model, which allows us to identify the many-body 
expressions for collective inertia B, deformation force /, and friction coefficient 7. 

Our microscopic expression M, ( !56l) . for inertia B has the form of inertia of collective 
motion in an abstract fluid with the velocity Ux = 5'V0x = '^"^Vxx- It should be stressed that 
—m~^Xx has nothing to do with the potential of the hydrodynamic velocity simply because this 
latter does not make sense in the nuclear case, once the nucleon mean free path A is comparable 
or greater than the nuclear sizes. To the best of our knowledge this work is the first to derive 
the expression f = Ag for the deformation force, where is the Lagrange multiplier used 
in the expression for pg to constrain it to a given expectation value q of collective coordinate 
Q. Having shown that A^ possesses the property A^ = —dgF, where F is the free energy, we 
establish another expression for the deformation force: / = —dqF, which though used in some 
analyses of collective motion, has never been derived from first principles. Our microscopic 
expression for 7, fll26p . relating this quantity with the time integral of a time correlation 
function of (p = J (ix^x^x for the pg ensemble, is also novel. 

For practical calculations of /, M, and 7 the real nucleus is replaced with a self-confined 
Fermi gas (SCFG) of quasiparticles, which, unlike the standard gases confined by external 
forces, is confined by its own forces. The SCFG model is based on the postulates that the one- 
body density matrix and the entropy of strongly-interacting nucleons, represented by Pg have 
the same form as those of noninteracting fermions bound by an external one-body potential 
and that the expectation of the bare interparticle interaction for pg may be replaced with 
the expectation of the effective interaction for the two-body density matrix of noninteracting 
fermions. Use of the extreme property of Gibbs grand potential constructed within the SCFG 
model for the pg ensemble allows one to derive the temperature-dependent constrained Hartree- 
Fock equations with the effective interaction. These equations become especially simple when 
the effective interaction is the Skyrme force. 

To find a SCFG expression for 7 from Eq. fll26p . we assume independent decay of particle 
and hole states. Replacing for simplicity the decay width Fa with the effective width F not 
depending on quantum numbers, we arrive at the formula 7 = 2vpM/ A, where Vp is the Fermi 
velocity, M is the collective inertia, and X = Vp/T the mean free path of the nucleon. Accord- 
ing to this formula the A, q, and temperature dependences of the ratio 7/M are completely 
determined by those of fp/A = F. Perhaps this result will foster further analytical and com- 
putational studies of the nucleon mean free path A in spherical, deformed, and superdeformed 
nuclei for realistic versions of the Skyrme force. 

One of the principal results of this present paper and our previous work [35] is that they 
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stress that the hnear response theory (LRT) expressions for B and 7 [27l[28l[29| [30 t[3T|[32l[33l [M] 
are wrong. This is because the LRT is not apphcable in the nuclear case. Namely, one cannot 
treat the force distorting the shape of a time evolving nucleus as an external force, because this 
force is depending on the intrinsic state of the nucleus. Moreover, the density matrix of the 
LRT implies that the mass flow in the bulk of the time-evolving nucleus is absent. 

To understand the nature of nuclear collective motion it was important to realize a key role 
of the mass flow in this phenomenon and obtain an explicit expression for this quantity in a 
slowly evolving nucleus: if the equilibrium state of the nucleus is represented by the real-valued 
single-particle wave functions ipx (x), for which the mass flow 

h = - (V2) (<^IV^A - (V<^^) ^x) = 0, 

then in the quasiequilibrium state of the nucleus those wave functions acquire the form e^^'^ipx (x), 
to give rise to a nonvanishing mass flow 

jx = XI f'^'^l W ^^'^ 

A 

where fx are the Fermi gas occupation numbers. 

It is interesting to note that this way of generating the collective motion in hot nuclei 
is similar to generating the coherent motion in the entrance channel of the head-on nucleus- 
nucleus collision leading to the composite system under study (the total angular momentum 
of the nucleus is in our model). Really, the ground states of two colliding nuclei can always 
be represented by superpositions of Slater determinants of single-particle standing waves. The 
relative motion with momentum in the head-on Ai + Ai collision along the X3 axis then gives 
rise to the phase factor e^"^^ attached to those standing waves with 

V?x = -{PrjA^xj, if X3 < X3 and v^x = (PrM2)a;3 if X3 > X3, 

where X3 is the touching point. This new vision of mechanism of collective motion as manifes- 
tation of coherency of the phases of individual nucleons in the hot nucleus qualitatively differs 
from the N. Bohr picture of collective motion in flssion [Ij, according to which the collective 
motion arises due to coherency of the velocities of individual nucleons. 
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Appendix A 



In this appendix the operators Q, Ai, P, K, H, Q, Ai, P, He, V, K'^ are expressed in terms 
of the appropriate primed operators. 

Using nx = n'^ in ( !58l) . ( !60|) . we obtain 



Q = Q', 
M = M'. 

Eq. (!59l) for P on accounting for (|28ll . ( I48l) takes the form 

P = P' + qM', 

where 

P' = J (ixUxPx, = m y dxu^n'^. 

Eq. (l89i) for after accounting for (11671) . ( 11661) takes the form 

K = P'. 



Eq. ( 1521) for 7i on accounting for (H8l) . (I168p becomes 

n = n'+ qP'+^q^M'. 



Employing Eqs. (ITTOll . (1T651) in Q = i [H, Q], we use [M', Q'] = to obtain 

1 



Q = i 



n'+ qP'+^q^M',Q' 



I [w + qp', Q'] = g' - gir;, 



where K[ = i [Q', P']. The equality [Ai', Q'] = follows from the relation J dyfy 
valid for an arbitrary function fy. Employing Eqs. (I170|) . (I166|) in At = i [H, Ai], we find 



(165) 
(166) 

(167) 
(168) 

(169) 
(170) 

(171) 
0, 



M = i 



n'+ qP'+-q^M',M' 



M' - qK'^ 



2? 



where K'^ = i [At', P']. Employing Eqs. (IT70|) . (11671) in P = z [H, P], we have 

1 



P = i 



n'+ qP'+^q^M',P' + qM' 



(172) 



(173) 



Eq. ([70]) for on accounting for Eqs. (IT70D . (IT671) . (IT65D . (IT66D takes the form He = n'+ 
AgQ', or using Eq. ([71]) for Hg, 

He = H'g. (174) 

Eq. (I7T1) for on accounting for (11651) . (11671) . (I166p becomes 

1 



V = AgQ' - qP' - ^q'M'. 

2 



Eq. (imp for on accounting for (11740 . (I169p takes the form 

K' = i[n'g,p'] =T', 

where T' is the primed partner of T , introduced in 



(175) 



(176) 



33 



Appendix B 



In this appendix we derive the integral equation fl971) for Kt 
derivation is taken from [42j. Define the operator 



where = e^^^-ft'e On using; the identities 



J{t) = J(0) + [J{t) - J(0)] = J(0) + / dsdsJ{s), 

Jo 

j(o) = k, 



9, J(s) = 9,e-*"''^e*""ire -^e 
where V = He — Ti., Eq. (11771) becomes 



JsHc 



(177) 



K+ / rfse 



sHc 



which leads in an evident way to the desired integral equation for Kt 



178) 



Appendix C 

Using the relation dtkn,t = (kr,,i^ for kn,t = {Kn,t)q and Eqs. §S>, ([79]), (02]), we find 

t„,, = /<;x«,„P.,,(x.,.). u>,,, = -u,Vr„ = -mu,VSj. (179) 

This equation says that kn,t are linear functions of t, which permits us to write dfkn^t = ^kn/r^, 
with 6kn = kn^T^ — knfl- Thus, taking into account that t] = r~^, we obtain 

/•oo . . roo poo 

/ dte-^Ukn,t) = / dte-'/^-dtkn,t= / dte-'/^-5kr,/T^ = 5kn. 
Jo \ Jo ^0 



From 6kn = kn,T„ — knfl and (11790 we find the expressions 



6ki = —6q ■ J (ixUxVyx<9gPx5 Sk2 = —Sq ■ m (ixUxVMx<9gPx^ 



;i80) 



where 5q = qr^ — q. Since 6q = qr^, Eqs. (I179p . (11801) show that dte '^^ (^Kn,t^ are 
proportional to q. 
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Appendix D 



Let us prove that {i [JF, Q]) g = 0. Exploiting the cyclic property of the trace and definitions of 
JF, Tig, and Q, we find that 

{z[J^,Q])=tT U \q,P 



Pg 

From definitions of Q, -P, Q and Eq. ( l44l) it follows that 



tr i 



Q,P 



Pg = tr U 



Pg 



-m-^ / c^xrfyF.V^M^tr (z pg) 



From the relation 



dxp dxa 

valid for any c-number function /y, and Eq. fll9l) . according to which tr(pxPg) 



0, it is seen 



that tr 



pg 



0, hence, {i [J^,Q])„ = 0. Analogously one shows that {i [J^,Ai]) = 0. 
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Figure 1: Equidensity lines for ^^°Yb calculated by the ETF method at the indicated values of 
the quadrupole moment (in barn). From outside to inside, and in units of po) the lines represent 
contours of constant density p = 0.1, 0.3, 0.5, 0.7, 0.9 and 1.1. The horizontal axis corresponds 
to the z direction and the vertical axis to the r direction (units in fm). 
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Figure 3: The volume of ^^'^Yb along the fission path in units of its volume at (^2 = from the 
ETF - SkM* calculation. 
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Figure 4: Comparison between semiclassical friction coefficients (dotted and daslied) and 
7 = 2f pM/ A (solid lines) along the fission path in ^°'^Pb. 
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Figure 5: Same as in the previous figure but for ^^^Ds. 
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